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Abstract. We present the results and parameterization 
of the 0.5-2 keV Luminosity Function of AGNs from var- 
ious ROSAT Surveys, ranging from the ROSAT bright 
Survey from the ROSAT All-Sky Survey (RASS) to the 
Ultra-Deep survey on the Lockman hole. A Luminosity- 
dependent density evolution model, where the density evo- 
lution rate drops at low luminosities, gives an excellent 
parametric description of the overall XLF covering wide 
ranges of redshift and luminosity. 

The number density evolution of high-luminosity 
AGNs in our sample shows a similar behavior to optical 
and radio surveys, except that we do not find evidence for 
the rapid decrease of the number density at z > 2.7. The 
discrepancy is marginally significant and including more 
deep survey results would make better determination of 
the behavior. 



1. Introduction 

Strong X-ray emission is a prominent key character of an 
AGN activity (we use the term "AGNs" for Seyfert nuclei 
and QSOs collectively). Thus unbiased X-ray surveys and 
identifications are important for investigating cosmologi- 
cal evolution of AGN activities. A combination of ROSAT 
surveys, ranging from the RASS to the ROSAT Ultra Deep 
Survey on the Lockman Hole, provides a large sample of 
soft X-ray selected AGNs. In this article, we report ba- 
sic results of our work on the soft (0.5-2 keV) X-ray lu- 
minosity function (SXLF) of AGNs and its evolutionary 
properties, using a combined ROSAT sample of about 670 
AGNs. 

A construction of a population synthesis model of the 
Cosmic X-ray Background (CXRB) with a combination 
of unabsorbed "type 1" and intrinsically absorbed "type 
2" AGNs are presented elsewhere (Miyaji et al. 1999; see 
also Schmidt et al. 1999). Unless otherwise noted, we use 
Ho = 50 h 50 [kms- 1 Mpc~ 2 ] and (fi m ,fi A ) = (1.0,0.0) in 
our calculations. The symbol L x refers to the luminosity 
in [erg s _1 ] in 0.5-2 keV using a K-correction assuming 
a power-law photon index of T — 2. This is equivalent 
to the no K-correction case, and more realisticaly, this 



should be considered as the 0.5(1 + z) - 2(1 + z) [keV] 
luminosity at the source, given the variety of AGN spectra. 
The symbol 5x14 refers to the 0.5-2 keV flux measured in 



[10 14 ergs 1 cm 



2. Sample 



We have constructed a combined sample from various 
ROSAT surveys as summarized in Table ||. RBS and SA-N 
are from RASS, RIXOS is a serendipitous survey and oth- 
ers are pointed deep surveys. For the UKD and Marano 
samples, we have included sources identified with QSOs 
and NELGs, but used the objects with 5 x i4 > 0.5, in 
order to minimize the possibility of including significant 
number of misidentified sources (see Schmidt et al. 1998). 
For the Lockman hole, we have combined the very deep 
(~ 1 Msec) HRI data (HRI offaxis < 12' with a limiting 
flux of = 0.19) and PSPC data (HRI offaxis > 12', 
= 0.38 - 0.5 depending on the PSPC off-axis angle). 
For the PSPC data, we have used only pulse-height chan- 
nels corresponding 0.5-2 keV and, for both the PSPC and 
HRI data, the countrate-to flux conversion has been made 
using a photon index of V = 2 corrected for the effect of 
the absorption in our galaxy. Unlike some previous works, 
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Fig. 1. The survey area of the combined sample are plot- 
ted as a function of the limiting 0.5-2 keV flux limit (solid 
line). For reference, [N(> S)} -1 for all the X-ray sources 
is overplotted (dashed-linc) in the same scale. 



2 



T. Miyaji et al.: ROSAT AGN Luminisity Function 



Table 1. ROSAT Surveys used in the analysis 



Survey a 


olim 
•J a:14 

[erg s" 1 cm -2 ] 


Area 
[deg 2 ] 


No. of 
AGNs 


RBS 


« 250 


« 20000 


223 


SA-N 


» 13 


w 640 


130 


RIXOS 


3.0 


« 15 


205 


NEP 


1.0 


« 0.2 


13 


Marano 


0.5 


w 0.2 


28 


UKD 


0.5 


w 0.16 


29 


RDS 


0.19-0.5 


« 0.1 


62 



a Abbreviations/Reference - RBS: The ROSAT Bright Survey 
(Schwope et al. 1998), SA-N: The Selected Area-North (Ap- 
penzeller et al. 1998, RIXOS: The ROSAT International X-ray 
Optical Survey (Mason et al. 1998), NEP: The North Ecliptic 
Pole (Bower et al. 1996) Marano: The Marano field (Zamorani 
et al. 1998), UKD: The UK Deep Survey (McHardy et al. 1998), 
RDS: The ROSATDeep Survey on the Lockman hole (Hasinger 
et al. 1998; Schmidt et al. 1998) 

we have included all objects which have AGN indications 
(except BL-Lac objects), not only the ones with appar- 
ent broad lines. Some of these objects would have been 
classified as "Narrow Emission Line Galaxies" (NELGs) 
by other groups and thus would have been excluded from 
their analysis. 

Fig. [j] shows the available area of the combined sample 
as a function of the limiting flux with the inverse of the 
AGN N(> S) in our sample. 



The LDDE expression we have used for the overall 
SXLF is: 

d / T (Z/X : Z) = A [(£ X /Z*) 71 + (W^PF 1 • e(z,L x ), (1) 
a .Log L x 

where e(z, L x ) is the density evolution factor: 

C (1 + z )max(0,pl-a(Log L a -Log L x )) ^ < ^ ^ < £j 

< (l + z)P 1 (z < z c ; L x > L a ) . 

[e(z c ,L x )[{l + z)/(l + z c )] p2 (z>z c ) 

The parameter a represents the degree of luminosity 
dependence on the density evolution rate for L x < L a . 
The PDE case is a = and a greater value indicates 
lower evolution rates at low luminosities. 

The best- fit parameters/90% errors for the LDDE 
model, 2D K-S probabilities P2DKS (probabilities that 
the D-value for the 2D K-S statistics exceeds the ob- 
served value) are shown in Table |[ The integrated 
0.5-2 keV intensity of the model, extrapolated below 
the survey limit, is also shown. This can be com- 
pared with the extragaxlactic CXRB intensity of (7.4- 
9.0)xl0 _12 [ergs _1 cm -2 deg -2 ], from an update of the 
ROSAT/ASCA measurements by Miyaji et al. (1998a). 

For demonstrations of the goodness of the overall fit, 
we compare the "flattened" Log N - Log S (S 1 - 5 N(> S)) 
and cumulative /(< z) curves of the real data with the 
model prediction in Figs. ^|and^ respectively. The results 
of the 2DKS test, along with these comparisons show that 
our overall description is an excellent representation of 



3. Description of the overall SXLF 

Using the maximum-likelihood method, we have fitted the 
unbinned sample with a number of SXLF models, includ- 
ing the Pure Luminosity Evolution (PLE), Pure Density 
Evolution (PDE), and the Luminosity-dependent density 
evolution (LDDE) models. The overall fit has been made 
for the redshift range 0.015 < z < 5 and 41.7 < Log L x < 
47.0. We have tested the best-fit models in each class with 
a two dimensional Kologomorov-Smirnov (2D K-S) test 
(Fasano & Franceschini 1987). The details of the statis- 
tical analysis will be discussed elsewhere (Miyaji et al. 
1998b). The basic results are: 

— Unlike previous works (e.g. Boyle et al. 1994; Page 
et al. 1996; Jones et al. 1997), PLE was absolutely re- 
jected. The 2D K-S probabilities for the best-fit PLE 
models are less than 10~ 2 for all sets of cosmological 
parameters considered (see Table. ||). 

— The best-fit PDE was marginally rejected with a 2D 
K-S probability of « 5%). One difficulty of the best-fit 
PDE model is that it overproduces the extragalactic 
soft CXRB intensity. 

— The LDDE model, where the evolution rate drops at 
lower luminosities, gave an acceptable overall descrip- 
tion of the SXLF (see below) . 
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Fig. 2. The l/V a estimates of the SXLF in different red- 
shift bins (as labeled) are shown. One sigma error bars 
are shown for each bin. Important upper limits are shown 
as downward arrows. The top tick and symbols mark of 
the upper limit correspond to 2.3 objects (90%) and one 
object per bin respectively. The solid lines show rough pre- 
dictions from the best-fit LDDE model (see below) for the 
redshift bins. 



Table 2. The best- fit LDDE Parameters 



(fim,fiA) Parameters/2DKS probability /Intensity 
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Flattened AGN N(>S) Curve (0.5-2 keV) 



(1.0,0.0) A = { 



1.57 ± .11) x 10 -e ; L* = 0.57+ g 
7i = 0.68 ± .18; 72 = 2.26 ± .95; pi = 5.4 ± .4 
2 C = 1.51 ± .15; p2 = 0.0 (fixed) 
a = 2.3 ± .8; LogL a = 44.2 (fixed) 
P2dks=51%; /xi2 = 4.5 ± 0.7 



(0.3,0.0) 



A 
71 



0. 



0+.21 

°-.13 

pi = 5.8 ± .4 



(0.3,0.7) 



;2.18±.15) x 10 -6 ;L« = 
0.58 ±.20; 72 = 2.26± 
z c = 1.54 ± .15; p2 = 0.0 (fixed) 
a = 2.6 ± .6; LogL a = 44.5 (fixed) 
P2dks=34%; 7 xl2 = 4.8 ± 0.7 



A = (2.04 ± .14) x 10 -6 ; L, = 0.50i ?g 
71 = 0.58 ± .19; 72 = 2.27 ± .08; pi = 5. 
z c = 1.51 ± .14; p2 = 0.0 (fixed) 
a = 2.9 ± .6; LogL a = 44.6 (fixed) 
P2dks=38%; Jxi2 = 4.6 ± 0.7 



± .4 



(0.0,0.0) 



-«;£.= 0.45t:S 



A = (2.35 ± .16) x 10" 
71 = 0.56 ± .20; 72 = 2.24 ± .08; pi = 5.9 ± .4 
z c = 1.55 ± .14; p2 = 0.0 (fixed) 
a = 2.7 ± .6; LogL a = 44.6 (fixed) 
P 2 dks=26%; 7xi2 = 5.3 ± 0.9 



Units 

[10- 12 



A: 

ergs - 



[hlo 



Mpc~ 
! deg" 2 



3 ], W. [10 44 
in 0.5-2 keV 



/i 50 2 erg s" 



the current data. The best-fit LDDE models, when ex- 
trapolated below the survey limit, gives 50-70% of the 
extragalactic background in the 0.5-2 keV band consider- 
ing errors of the fits. If the LDDE parameters are adjusted 
to produce 90% of the extragalactic 0.5-2 keV background 
(clusters should contribute ~ 10%), the 2D KS probability 
drops to ~ 8%. 

4. Evolution of QSO Number Density 

Number density evolution of luminous QSOs is one of 
the important pieces of key observational information on 
understanding blackhole formation and accretion history 
In particular, we put emphasis on high luminosity QSOs 
(Log L > 44.5), where the XLF is consistent with the 
slope of 7 = 2.3 at all redshift and free from complicated 
luminosity dependence of the evolution rate and contam- 
ination from star formation activity. 

The comoving number density of the luminous QSOs 
in our sample are plotted as a function of redshift in Fig. 
|. For comparison, number densities (normalized to be ap- 
proximately equal to the soft X-ray point at z ~ 2.5) of 
optically-selected (Schmidt, Schneider, Gunn 1995, here- 
after SSG95) and radio-selected (Shaver et al. 1997) QSOs 
plotted as a function of redshift (see caption). 

Unlike the optical and radio cases, we do not find evi- 
dence for decrease of the space density beyond z w 2.7. Us- 
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Fig. 3. The S 1 ^N(> S) (a horizontal line corresponds 
to the Eucleadean slope) curve for our sample AGNs is 
plotted with 90% errors and are compared with the best- 
fit PLE (short-dashed), PDE (dot-dashed), and LDDE 
(dashed) models. The dotted fish is from the fluctuation 
analysis of the Lockman hole HRI data (including non- 
AGNs) by Hasinger et al. (1998) 
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Fig. 4. The cumulative intensity /(< z) of the sample 
AGNs, defined by J2 Zt <z S X i/A(S xl ), where 5 X i is the flux 
of the object i and A(S X i) is the available survey area at 
this flux (Fig. is plotted as a function of redshift (thick 
solid line with 90% errors). The best-fit LDDE model is 
overplotted for the portion represented by the sample (S > 
S\ im = 1-9 x 10~ 15 [ergs -1 cm -2 ]; dashed line) and the 
total including extrapolation to lower fluxes (thin solid 
line; with a 90% fitting error). Also the range of the 0.5-2 
keV extragalactic background intensity by Miyaji et al. 
1998a (updated) is shown. 



ing the maximum-likelihood fitting, we have checked the 
significance of the inconsistency. Requiring that the num- 
ber density drops beyond z = 2.7 as SSG95, the likelihood 
value (varies as x 2 ) increased by 3.3, showing that the sig- 
nificance of the inconsistency is 93%, which is marginal. 
In the 3.0 < z < 4.6, we have 5 QSOs in the sample, while 
expected number in the presence of the decrease like the 
SSG95 result is 2.4. Including more surveys with a good 
completeness at the depth of S x ~ 5 x 10 -15 [erg s -1 cm -2 ] 
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• Soft X-ray {Log L„>44.5) 

A Optical. M B <-26.0, normalized (SSG95) 

* Radio, normalized (Shaver el al. 1997) 
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Fig. 5. The comoving number density of luminous 
(Log L > 44.5) QSOs in our ROSAT AGN sample are 
plotted as a function of redshift. For comparison, num- 
ber density of optically-selected (Mb < —26) (dashed line 
and filled triangles, from SSG95) and radio-selected (stars, 
Shaver et al. 1997) QSOs, normalized to the soft X-ray se- 
lected QSO number density at z ~ 2.5 are overplotted. Er- 
ror bars are at the 90% confidence level and upper-limits 
correspond to 1 object per bin at the symbol and 90% at 
the top. 



would enable us to discuss whether the soft X-ray selected 
QSO number density drops beyond z w 3. 

5. Discussion 

There are some discrepancies between our work and pre- 
vious ones from several authors from combinations of 
the Einsten Extended Medium Sensitivity Survey (EMSS) 
and ROSAT surveys (e.g. Boyle et al. 1994; Page et al. 
1996; Jones et al. 1997), in that our results are not consis- 
tent with PLE. This is partially due to the fact that our 
combined sample has more AGNs at fainter fluxes, thus 
better statistics excludes simpler description. However, a 
direct comparision of the 1/V a estimates of the XLF be- 
tween our sample and that of Jones et al., for example, 
shows no descrepancy in the z < 0.4 bin, and high lumi- 
nosity part (Log L x <; 44) of the z > 0.4, but descrep- 
ancics appear at the low luminosity end of the z > 0.4 
regime. The SXLFs in the regime of bright sources, i.e. 
the RBS/SA-N based on the RASS for our work and the 
EMSS sample used by them, arc mutually consistent. Thus 
the uncertainties in conversion of fluxes between ROSAT 
and Einstein bandpathes do not make significant contribu- 
tions to the descrepancy. Two likely major causes are: (1) 
inclusion of apparently narrow-line objects which have in- 
dications of AGN activities in our work, while they mainly 
argue the XLF for broad-line objects. (2): because of the 
misidentification of the faintest X-ray sources (Sxi4 Si 0.5) 
to field galaxies in their PSPC positioning, they miss some 
AGNs. A more complete analysis (e.g. comparing redshift 



distribution) using the original catalogs will be discussed 
in Miyaji et al. 1998b. 

The extrapolations of our best-fit LDDE models ex- 
plain 50-70% of the 0.5-2 keV extragalactic CXRB, con- 
sidering errors of the fits. Clusters of galaxies are expected 
to contribute about 10% in this band. The remaining con- 
tributors could be some low-luminosity galaxies/ AGNs 
(Log L x ^ 10 41 - 7 ) (see Fig. 2a of Hasinger 1998). These 
have about the same local volume emissivity as the AGNs 
in our SXLF, but uncertainties up to by a factor of few 
may exist, since this low-luminosity galaxy XLF is based 
on a sample of galaxies nearer than 7.15 Mpc and the role 
of local overdensity can be important. Also practically no 
direct observational information exists for the evolution of 
these low luminosity sources. 

If the low-luminosity galaxies/ AGNs contribute signif- 
icantly to the remaining CXRB, these can be: (1): intrinsi- 
cally low luminosity AGNs (2) high-intermediate redshift 
obscured AGNs, which are expected to contribute much of 
the harder CXRB and redshifted into the ROSAT band. 
(3): star formation activity. 

It is also possible that the behavior of the AGN SXLF 
(Log L x 10 41 ' 7 ) at intermediate-high redshifts does 
not follow the simple extrapolation of the current LDDE 
model below the luminosities corresponding to the survey 
limit at that redshift. Especially our particular formula for 
the LDDE tends to make the N(< S) drop rapidly below 
the survey limit. Also the LDDE model- integrated inten- 
sity is sensitive to the lowest flux objects in the sample, 
where incompleteness may play a role. These are funda- 
mental uncertainties in the integrated intensity extrapo- 
lated from the best-fit representations of some functional 
form. 
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